Quantum photon emission from a moving mirror in the nonperturbative regime 
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We consider the coupling of the electromagnetic vacuum field with an oscillating perfectly- 
reflecting mirror in the nonrelativistic approximation. As a consequence of the frequency modulation 
associated to the motion of the mirror, low frequency photons are generated. We calculate the 
photon emission rate by following a nonperturbative approach, in which the coupling between the 
field sidebands is taken into account. We show that the usual perturbation theory fails to account 
correctly for the contribution of TM-polarized vacuum fluctuations that propagate along directions 
nearly parallel to the plane surface of the mirror. As a result of the modification of the field 
t-O , eigenfunctions, the resonance frequency for photon emission is shifted from its unperturbed value. 

<D . 

■ PACS: 4250, 0365, 4210. 

C\l ' I. INTRODUCTION 

> : 

When the boundaries of the quantum radiation field are set in motion, low frequency photons may be excited out 
of the vacuum field state. Such quantum effect is closely related to the radiation pressure exerted on the moving 
_^ ' boundaries by the vacuum field jjj. Given its conceptual importance, one has proposed an experiment to observe 
the effect with the aid of high-Q microwave cavities J2|. Exact relativistic results are known for one-dimensional 
OS ' models Q, but the analysis of realistic 3D configurations usually demands some perturbation method. Several 
different perturbation approaches Q- [Q have been successfully employed to compute the radiation effect induced by 
r"| • the nonrelativistic motion of a broad variety of 3D boundaries. In reference ||, the photon emission rates induced by 
O || the oscillation of a single perfectly-reflecting plane mirror were obtained by taking the long wavelength approximation 
, 1 1 . and assuming the motion induced corrections to be small. More recently, this particular perturbation approach was 
also employed in the problem of two plane mirrors forming a cavity |7) . 

The connection between the nonrelativistic limit and the long wavelength approximation may be understood by 
considering the example of a mirror oscillating at frequency Slo with amplitude Sqo. Then the nonrelativistic limit 
translates into Qofiqo ^ c, and since the emitted photon frequencies satisfy the inequality uj < Qo ®: they correspond 
to wavelengths A much larger than Sqo. In this long wavelength regime, the scattering of a plane wave of frequency 
uj usually generates sidebands at frequencies ui — mQo with amplitudes of the order of (S/o/A)'™', where m is any 
integer. Hence the sidebands have amplitudes which are increasingly small as the order m increases. Moreover, the 
motion induced correction of the elastic scattered component is usually of the order of (Sqo/X) 2 . Therefore, the photon 
emission effect induced by a nonrelativistic motion may be analyzed by taking a first-order expansion in Sqo/X, and 
then the entire motion effect is contained in the sidebands at w±Sl . Accordingly, in this approach the motion induced 
terms are treated as small perturbations which may be computed directly from the field for a mirror at rest 

However, when the incident wave is polarized so that the magnetic field is perpendicular to the plane of incidence 
(TM polarization), the first downshifted sideband (m = —1) may be intense if the parameters are such that it 
corresponds to a grazing wave (peak amplitude value of the order of X/Sqo). Such enhancement effect, which is 
reminiscent of the Wood's anomalies in the scattering by a diffraction grating at rest, first explained by Rayleigh in 
the beginning of the century jq], was analyzed in reference ^] from the point-of-view of classical wave optics, always 
in the long wavelength approximation. It takes place over a very narrow spectral interval, and may be interpreted in 
terms of resonances associated to the plane symmetry ( in the context of diffraction gratings, a discussion along these 
lines was proposed by Fano (l(]]). More specifically, the resonances are related to surface waves (plasmons) which 
are the field eigenfunctions for the unperturbed system (surface at rest). In the perfect-reflection limit considered 
here, they are expelled from the medium and then degenerate into grazing travelling waves. These unperturbed 
eigenfunctions are TM polarized (so that the electric field is perpendicular to the surface of the mirror), thus no 
resonance enhancement occurs for TE scattering in the long wavelength approximation. 

As for the motion induced quantum radiation effect, in order to compute the emission rate taking these resonances 
into account, we follow in this article a nonperturbative scheme, though still assuming a nonrelativistic motion, which 
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allows us to employ the long wavelength approximation. We consider in detail the contribution of TM-polarized 
vacuum fluctuations which propagate along directions nearly parallel to the surface of the mirror. As we show below, 
the motion induced correction to TM-polarized nearly grazing waves is not a small perturbation in this case, and the 
dynamical coupling between different sidebands must be taken into account. As a result the elastic field components 
are strongly modified, leading to a shift of the resonance frequency for generation of photons. 

The paper is organized as follows. In section 2, we expand the scattered field as a superposition of sidebands in 
order to derive a nonperturbative solution of the boundary condition associated to an oscillating mirror, which is then 
used in section 3 to compute the photon emission rate. Section 4 presents a discussion and some final remarks. 



II. BOUNDARY CONDITIONS 



We consider a plane perfectly-reflecting mirror oscillating along the direction perpendicular to its surface, which 
we take along the x axis. The mirror's position at time t is given by the equation x = Sq(t). We use MKS units 
and take eo = 1, c = 1. As usual we decompose the electromagnetic field into components according to whether 
the electric field is perpendicular (TE) or parallel (TM) to the plane of incidence. As a consequence of the plane 
symmetry, the two polarizations are uncoupled in the scattering by the moving mirror, and accordingly we have two 
separate problems to solve. In the case of TE polarization, the perturbation approach is always valid in the long 
wavelength regime (as in the problem of reflection by a shallow diffraction grating) , so that we only consider in this 
paper the generation of TM photons out of the vacuum radiation field. We define the vector potential A as follows: 

e (tm) = v x ^ (1 ) 

B<™> = d t A. (2) 

Moreover, it satisfies the Gauge condition 

V-.A=0. (3) 

The key motivation for defining the vector potential as in eqs. is the fact that A ■ x — 0, which entails that 

A does not change under a Lorentz boost. For the Lorentz frame in which the mirror is instantaneously at rest the 
boundary condition for A simply states that its normal derivative vanishes at the mirror's position. Then, translating 
back to the laboratory frame, such condition yields (see for details): 

[d x + Sqd t ]A(Sq(t),T h t)=0, (4) 

where r|| = (y, z). 

The field is written as 



A = An + &Arct (5) 

where &/4 rc t is the retarded field that represents the effect of the motion of the mirror, whereas Ai n , representing the 
input field (t — > -co), is itself a sum of incident and reflected waves, and hence satisfies the boundary conditions of a 
mirror at rest. We take periodic boundary conditions on the plane of the mirror over a square of surface S, which is 
to be identified with the surface of the (very large) mirror later. Therefore, the normal mode decomposition of A\ n 
is given by 

A iu (x, r,| ,t) = J o °°^-J2 y^cos^zK 1 * 11 "" 1 ' e- ifc "*a inn (A: x )e„ + H.c, (6) 

where H. c. means the Hermitian conjugate, ky = (n y y + n z z)2"K/\~S, with n = (n y , n z ) denoting a pair of integer 

numbers, e„ = x x k|| n is the polarization vector, and k n = \Jk% + fcjj^ is the frequency of a given normal mode. The 
input bosonic operators satisfy the commutation relations 

[a inn (k x ),a ian ,(k' x )] = (7) 

and 

[ainn(kx),amn'(K)^} = ^5(k x - k' x )S n , n : . (8) 
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We assume that the oscillatory motion of the mirror is nonrelativistic: Sq <C 1. As discussed in the introduction, 
this assumption allows us to take the long wavelength approximation in order to solve for &4. re t in terms of *4.i n - Thus, 
we expand the fields around x = keeping terms up to order of Sq/X, and find, from eqs. (||) and (0), 

(8 X + 8qd 2 x + Sqdt) SAret (0, r„ ,t) = - (5qd 2 x + Sqdt) An(0, r,| , t) (9) 

In the perturbative regime considered in reference ||, SA. rc t is already of first-order in Sq/X, and then only the first 
term in the lhs of eq. (||) contributes to the required order of Sq/X. In this paper, we consider a regime, involving 
propagation near the grazing direction, where <5A re t may be of a higher order over a very short spectral range. In 
fact, eq. (|^) itself provides a hint that it is not possible to neglect the last two terms in its lhs when considering waves 
propagating along the grazing direction, since d x SA rc t vanishes in this case. These terms provide a coupling between 
the motion induced sidebands, which is to be analyzed later in this section. 

It is convenient to define a Fourier representation that takes advantage of the plane symmetry and the corresponding 
conservation of the parallel component of the wavevector, k|| . Hence we define the Fourier components of a vector 
field F(r,t) as 



(10) 



F n [x,u] = ^ J dt J s d \ z lUt e- 4k ""- r "F(x,r||,i), 
and the wave equation reads 

d 2 x F n [xM = (h 2 n -^Wn[xM- (ii) 

The normal mode expansion of eq. (ra) is written as 



where 



Ainn[xM = 6(/C(o;) 2 )W ^j^L cos(ACM:e) [e(w)o in „ - 0(-w)4,_ n (-£("))] e n , (12) 



K(u) = lim [{io + ief - fcjf ] 1/2 (13) 

e— >0+ 



is a function in the complex plane of u) with a branch cut along the interval in the real axis between — fcii and fcii. It 
represents the x-component of the reflected wave associated to the parameters u> and ky . The factor 0(/C(w) 2 ) in 
eq. ( [l2] ) (where Q denotes the Heaviside step function) is a consequence of excluding evanescent waves from the mode 
expansion of eq. (||), whereas the factor 0(w) establishes the connection between annihilation (creation) operators 
and positive (negative) frequencies, which is crucial for the understanding of the photon emission process discussed 
in this paper. 

The mirror oscillates around the position x = : 

5q(t) = (SftCOslV, (14) 

with f2o<%) *C 1, leading to the generation of sidebands at frequencies 

LO m = lu — mflQ. 

We take the Fourier transform of eq. (JoJ) as defined by eq. (|l0|) and use eq. ( |TT| ) to find (from now on we omit explicit 
reference to k| ): 

d x SA Ict [0,uj} + ~6q [( k \\l - (&M0,Wj] + A„[0,^])] = 0. (15) 

i=-i,i 

In order to solve eq. ([l5j), the retarded field is written as a superposition of all sidebands generated from the input 
field. In the half-space corresponding to x > 0, we have 

00 

SA ret [x,Uj} = exp (HC(lu)x) ^ 9m{U-m) -Aiu [0, OJ-m] ■ (16) 



m— — oo 
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According to eqs. ([L3|) and (|16|) (and to the particular choice of the branch cut on the complex plane associated with 
the former), 5A rc t[a;,a;] either corresponds to a wave propagating from the mirror into the half-space x > (when 
| a; | > k) or to an evanescent wave decaying along the positive x direction, since K, is of the form K.{pS) = i\K(u))\ when 
M < k. 

We replace eq. ([l6]) into eq. ([l5]) and use eq. ( |l3| ) to find, after changing the summation index: 
°° f <fyo 

^2 \ ifc(u)g m (u>- m ) — ^- {[K.{uj)) 2 + jujjD, )(g rn+j (uj- m } + 6 m ,-j) } A in [0,uj- m } = 0. (17) 

m— — 00 I j= — 1,1 



Since eq. fllTj ) applies to any input field, the expression within curly brackets must vanish for all values of m and 
arbitrary values of u>. When considering a given value of m in eq. ([l^), we replace u by u> m — u> — muio in order to 
have the g functions evaluated at the same frequency uj. We obtain an infinite system of coupled linear equations for 
the functions g m (oS) of the form 

Mg = Y, (18) 
where M is the symmetric tridiagonal infinite matrix given by M m m = iIC(u) m ), 

M m+ i jin = M„ V(Tl+ i = —H(u> m ), 
and with M m ,„/ = otherwise. We have introduced the auxiliary function 



H{uj) = loQq — K,(oj) 2 , 



and Y and g are column vectors, with 



while the elements of g are the functions g m (oj), the integer index m running from — 00 to 00. We may recover the 
perturbation results by neglecting the nondiagonal elements of the matrix M, thus resulting in uncoupled equations 
whose solution is 

gi{uj) ra 2 ^- s - l(u) M 2 ^l^JTy (19) 

and with all other functions, including Qq{lo), of higher order of LoSqg. Note that these results are in agreement with 
the analogous expressions of reference Pj , and could be more easily obtained directly from eq. (||) by assuming 5A rc t 
to be a small perturbation. 

We seek a solution that interpolates this perturbative regime, considered in detail in reference || , with the regime 
corresponding to the scattering of vacuum fluctuations propagating along the grazing direction, always in the long 
wavelength approximation. Thus, we assume that the x component of the input field is small: |/C(w)| <C u 2 8qo- In this 
case, it is no longer possible to neglect the nondiagonal terms in the row corresponding to m = 0, so that we cannot 
neglect the coupling between go, g\ and g-\. We then solve eq. (|l8|) with m running from —1 to 1 



qi(ui) — ——, — r o , (201 



and 



/C(oji) K(uj— 1) 

/C(o;i)g(tj-i) 

g - i(a,)= ]Q^)m 9li) (21) 

JC(oj-i ) 

g (u) + l = -2i ^/,v giH. (22) 

The perturbative results as given by eq. ( |l9| ) are recovered from eqs. (pp|)-(p^) when |/C(o;)/cj| ^> (cj<5go) 2 - On the 
other hand, since go is of the order of one when \tC(u>)/u)\ < (ujfiqo) 2 , the oscillation of the mirror generates a strong 
modification of the elastic field components when considering propagation close to a grazing direction. In fact, eq. ( [22] ) 
yields go{u) = —1 when K.(u>) — 0. Thus, in this limit, we have A. = An + &Aret — 0, showing that in this case the 
retarded field &A rc t is of course not a small perturbation |l2), and that the TM grazing waves, which are the field 
eigenf unctions in the case of a perfectly-reflecting mirror at rest, are no longer allowed in the motional case. This 
effect is at the origin of the resonance frequency shift for quantum photon generation out of the vacuum field state, 
to be discussed in the next two sections. 
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III. PHOTON EMISSION RATE 



In this section, we use the nonperturbative results for the retarded field found in section 2 in order to consider the 
effect of photon generation induced by the motion of the mirror. As discussed in connection with eq. ([l2]), changing 
from positive to negative frequencies yields a coupling between annihilation and creation photon operators, which is 
responsible for the radiation effect considered in this paper. 

In order to analyze such coupling, we start by writing the total field, as given by eq. (^), in terms of the advanced 
solution 5/4ad v °f the boundary condition. In this case, the homogeneous component of the solution is interpreted as 
an output field, which represents the limit t — > oo : 

Ai n [x,oj] + §A iet [x,u>] =A out [x,u)]+&A a ,dv[x,u]. (23) 

The output field .4o U t may be expanded in terms of output bosonic operators (also satisfying the commutation relations 
of eqs. (0) and (||)) and normal modes exactly as in eqs. @ and ([l^), because it also satisfies the homogeneous boundary 
condition associated to a mirror at rest. Moreover, the advanced solution may be expanded as in eq. (fl6|): 

+oo 

<i4. a dv[a;,a;] = exp (-iJC(u})*x) ^ fm(^-m)-Amt[0, ui- m ], (24) 

m— — oc 

where the coefficients f m are given by f m — g* n . 

Eq. (|23] ) establishes a connection between input and output bosonic operators. In what concerns the dependence 
on the spatial coordinate x, note that the evanescent components contained in both 6A. TC t and &/4 a dv cancel each other 
in eq. ( P3|) and do not contribute to the photon emission effect, as expected since the radiation field is associated 
to travelling waves only. Then, we take \uj\ > ku n and assume K.(u>) to be real in what follows. As a consequence, 
the fields may be written entirely in terms of the two linearly independent functions exp(i)C(uj)x) and exp(— ilC(uj)x). 
Since eq. ( |23| ) applies to any value of x, the coefficient multiplying exp(ifC(cj)x) must vanish, yielding, with the aid of 
eq. ©: 



-0(-w) [a out _„(-/C(w))-a in _„(-X;(w))] t }e n = 44 ret [0,a;]. (25) 



If we take to < —ku n we single out the creation operators in the lhs of eq. (|25|). However, the retarded field in 
its rhs evaluated at u> may also contain annihilation operators, since u> may correspond to a downshifted sideband 
associated to an initially positive frequency lo + f2o- Such mixture between creation and annihilation operators entails 
that photons are created from vacuum by means of the frequency modulation associated to the motion of the mirror. 
In Appendix A, we start from eq. ( |25] ) and compute the average value of the output number operator associated to 
given values of k|| and k x over the input vacuum field state |0 in) (we use the shorthand (...) to denote such average): 

( a t utti(A;x )a outn (M) = ^e(fio-fc- fciij '^oo-fcf ^' (26) 



where At is a coarse-grained time interval, and k = + k^ is the photon frequency. Eq. ( J26| ) explicitly relates the 

photon creation effect to frequency downshifting, which is dominated by g\ in the long wavelength approximation (an 
equivalent description in terms of frequency upshifting is also possible and leads to the same final results). The step 
function in its r.-h.-s. expresses the requirement that the input vacuum fluctuations at frequency flo — k correspond 
to traveling waves. We consider a direction of emission forming an angle 9 with the x axis, so that kn = fcsin#. 
We analyze the behaviour of the photon emission rate at fixed frequency k and propagation angle 8 as we tune the 
mechanical frequency Qq. It is convenient to define the dimensionless variable 

A^-l-sin#. 
k 

According to eq. (|26|), no photons are generated when A < 0. At such low values of the mechanical frequency, there are 
no travelling wave vacuum field modes at frequency f^o — k that match the requirement associated to the conservation 
of the parallel component of the wavevector, kn . As a consequence, regardless of the choice of 9, the mechanical 
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frequency must be higher than the analyzed photon frequency k. As could be expected, high frequency field modes 
are unaffected by the motion of the mirror (quasistatic limit). This basic fact is essential to the long wavelength 
approximation employed in the paper. 

At A = 0, the vacuum fluctuations that contribute to the photon emission effect correspond to grazing waves. 
More generally, when A <C 1 the wavevectors of the relevant vacuum fluctuations have their x component given by 



/C(O -k) = kV2 sin 9 A [f + 0(A)] . (27) 

In the perturbative approximation, the function gi remains finite (see eq. (U9h) and then the photon number as given 
by eq. (^6|) diverges as A — > since JC(Hq — k) vanishes in this limit. Ref. || provides a simple interpretation of the 
divergence at A = 0. In the perturbative theory, photons are emitted in pairs. When the mechanical frequency is set 
to A = 0, the observed photon is such that its 'twin' propagates along a grazing direction Jl3| . TM polarized grazing 
waves are eigenfunctions for the perfectly-reflecting mirror at rest — in fact they represent the limit of the surface 
plasmons of a conducting surface as we take the perfectly-reflecting limit. Thus, A = is the resonance frequency 
for the field given the boundary condition of a perfectly- reflecting plane mirror at rest, for in this case the external 
modulation perfectly matches the conditions for the creation of grazing photons. 

On the other hand, according to the nonperturbative result of eq. (p(i|), the coefficient g\(^lo — k) is proportional to 
/C(Hq — k) as A — > 0. Therefore, instead of the divergence at A = predicted by the perturbative theory, the photon 
emission rate vanishes in this limit. In Appendix B, we compute the angular distribution rate for A > from eq. (|2^): 

d 2 R S k 2 cos 2 9 , , , ,. . „. ll2 

1 01 [u = ft(A + sin0)]p. (28) 



dkdQ 2tt 3 ^A(A + 2 sin( 

It follows from eqs. (pp|), ( p7|) and ( |28| ) that the photon emission rate goes as \/A in the limit A <C (kSqa) 4 . On 
the other hand, according to the discussion following eqs. (p0|)-(|22]), the perturbative regime is recovered when the 
vacuum fluctuations propagate along a direction not close to a grazing direction: /C(f2o — k)/k 3> (kSqo) 2 , which from 
eq. ( p7[ ) translates into A ^> (kSqo) 4 . In this limit, we recover the results of Ref. 

Since g\ is dimensionless, it depends on k only through A and kSqo- Thus, apart from the trivial k 2 dependence 
in eq. (pq), we may replace Sqo, f2 arj d k by the two dimensionless parameters kSqo and A. In figure 1, we plot the 
photon emission rate as given by eq. (|2^) as function of A for 9 = 78° and kSqo = 0.03. In the main figure, we also 
plot the result from perturbation theory (dotted line), which is obtained from eqs. (|l^) and (|28|). As discussed above, 
whereas the latter diverges at A = 0, the exact emission rate (solid line) vanishes in this limit. However, the exact 
result also displays a singularity, which is shifted to a higher value A s . Thus, the resonance frequency for photon 
emission is shifted from its unperturbed value Qq = fc(l + sm9) by an amount 

<50 = fcA s . 

The exact value of A s is best displayed in the insert of figure 1, where we compare the analytical result as given by 
replacing eq. (^0|) into eq. (p8|) (dashed line) with the result obtained from the numerical evaluation (solid line) of the 
coupled equations ( Jl8| ) for g m with m running from —3 to 3 (we have checked that truncations at higher dimensions 
had not changed the results within machine accuracy). The insert shows that the exact frequency shift is slightly 
smaller than the analytical value for the particular values of 9 and kSqo considered in fig. 1. On the other hand, 
outside the neighbourhood of A s the analytical result works fairly well, and in fact for the scale employed in the main 
plot of fig. I the two methods provide indistinguishable curves. 

From eq. ( pp| ) , we may derive an analytical approximation for A s up to lower order of k&ja : 

A s = ^ sin 3 9(1 + sinfff { J- - 1 ) (29) 

32 V V3sin 2 6> + 4sin<2+ 1/ 

For 9 — 78°, and kSqo = 0.03, eq. ( |29| ) yields A s = 7.187 x 10 -6 , which is in agreement with the plot of the emission 
rate as given by the analytical result (dashed line in the insert of fig. I). In figure 2, we plot A^ as given by eq. ( |29| ) 
as a function of the observation angle 9. The shift vanishes at 9 = and increases by several orders of magnitude as 
9 increases from zero to 90°. In the next section, we discuss such behaviour starting from the analysis of the diagram 
representing the input vacuum fluctuations that contribute to the radiation emission effect. 
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IV. DISCUSSION 



In the previous section, we discussed the behaviour of the photon emission rate for fixed values of the photon 
frequency and direction of emission, as we tune the mechanical frequency f2o- In this section, we consider the com- 
plementary situation, in which Qq is fixed, in order to understand why the resonance frequency is shifted from its 
unperturbed value. For the sake of clarity, we analyze the photon emission process as an effect of frequency upshifting 
(rather than downshifting) in this section. 

Each TM-polarized photon have given values of frequency uj, parallel component of wavevector feii , and an azimuthal 
angle defining the direction of kn . The photons generated from a mechanical modulations at frequency FIq correspond 
to points in the ABC triangle (dark grey) in the uj x k\\ plane of figure 3, since they obey the inequalities uj + < Qo 
and uj > fcii > 0. Normal modes with uj < ku do not correspond to travelling waves and hence are not relevant here. 
They correspond to the light grey region in figure 3. 

In the uj x kit plane the frequency upshifting induced by the motion corresponds to a horizontal shift (since k\\ is 
conserved) by an amount of Slo- Hence the vacuum fluctuations that contribute to the emission effect correspond to 
the A'B'C triangle in figure 3. Vacuum fluctuations associated to the B'C line propagate along grazing directions. 
They give rise to photons associated to points along the BC line, all of them corresponding to A = 0. In this paper, 
we have shown that in the neighbourhood of B'C the perturbative approach fails to provide the amplitudes of the 
sidebands. Whereas the perturbative theory predicts a divergent emission rate all along the BC line, we have found 
that the emission rate vanishes in this case, and that the singularity is displaced into the interior neighbourhood of 
the BC border. 

This effect is a consequence of the motion induced modification of the field eigenfunctions. For a metallic surface 
at rest, the eigenfunctions correspond to surface plasmons. In the perfect-reflecting limit, they are expelled from the 
interior of the medium and degenerate into TM polarized grazing waves. In the framework of first-order perturbation 
theory, the dynamical modification of the field eigenfunctions is neglected, and then the resonant enhancement takes 
place when considering vacuum fluctuations along the B'C line, which give rise to photons in the field modes corre- 
sponding to A = 0. However, this resonant regime is inconsistent with the underlying assumption that the motion 
effect is a small perturbation. In this paper, we have employed a nonperturbative approach, and found the correction 
of the elastic field component, which is represented by the coefficient go in eq. (|l6|), to be very important in the case 
of nearly grazing waves (see eq. (|2^) and the discussion that follows). Thus, TM polarized grazing waves are no longer 
acceptable solutions in the motional case. Since the field eigenfunctions are strongly modified, the resonant region is 
displaced away from the B'C line in figure 3 into the interior of the A'B'C region. As a consequence, the critical 
region in the positive- frequency part of the diagram is displaced from the A = line (BC line), leading to a shift of 
the resonance frequency. 

Vacuum fluctuations associated to the A'B' border give rise to photons propagating along a grazing direction (AB 
line). In this case, the perturbative approach also fails to provide the correct photon emission rate (see remark |L2|]). 
The two critical regions overlap at point B, which corresponds to grazing photons at the subharmonic frequency, 
generated from grazing vacuum fluctuations. As we approach B along the BC line, the width of the nonperturbative 
region increases, as well as the frequency shift A s . In fact, that corresponds to the limit 8 — > 90° discussed in 
connection with eq. ( p9| ) (see fig. 2). Note, however, that eq. ( p9| ) was derived from eq. ( p0| ) by assuming A s <C 1. 
Moreover, eq. ( p0| ) itself was obtained from the coupled equations (18) by neglecting higher order sidebands, which 
turns out not to be a good approximation in this limit. As shown in the insert of fig. 1, the exact shift is smaller than 
the value given by eq. (|2^), and the discrepancy increases as 6 — > 90°. 

The ABC triangle is closed from below by the AC border. The latter corresponds to photons that propagate 
along the normal direction. Since they are generated from vacuum fluctuations that also propagate along the normal 
direction (A'C line), the nonperturbative theory provides no relevant correction in this case. As expected, according 
to fig. 2 the frequency shift A s vanishes as we approach the neighbourhood of point C (9 — > 0). 

In this paper, we have analyzed in detail the resonance effect associated to the excitation of surface plasmons 
from the vacuum field state when a mirror oscillates in free space. Due to the assumption of perfect reflectiveness, 
the surface plasmons degenerate into TM-polarized grazing waves, and the resonance appears as a singularity in the 
photon emission rate. The nonperturbative approach developed in this paper allowed us to calculate the shift of the 
resonance frequency induced by the motion of the mirror, which is a consequence of the dynamical modification of 
the field eigenfunctions. 

The authors thank Programa Especial de Visitante (PREVI-UFJF), Consclho Nacional de Desenvolvimcnto 
Cientffico e Tecnologico (CNPq) for partial support (J. P. R. F. M. and P. A. M. N.) and Programa de Nucleos 
de Excelencia (PRONEX), grant no. 4.1.96.08880.00 - 7035-1 (P. A. M. N.). 
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APPENDIX A: AVERAGE VALUES OF OUTPUT NUMBER OPERATORS 



In this appendix, we derive eq. (|26| ) for the average value of the output number operator. We replace the retarded 
field in eq. ( |25| ) by its expansion in terms of sidebands as given by eq. ([l6|) and use the normal mode decomposition 



as written in eq. (|12|). Then, as mentioned in sec. 3, we pick up a frequency to < — kn to find 



-K,(w))a out _ n (-K.(uj)) 









LU 




/C(w_i)_ 



e(w. 



\gi{u-i)y 



(Al) 



Note that it is also possible to compute the photon emission rate by picking up a positive to in eq. fl2q), but in this 
case the dominant contribution comes from the function <7_i, which describes frequency upshifting. Of course this 
method leads to the same final result for the average output number opera tor. 

In order to calculate the vacuum correlation function in the rhs of eq. ( |Al[ ), we introduce a coarse-grained time 
scale At, which satisfies QqM, ^> 1, by means of a frequency distribution sharply peaked around the mechanical 
frequency Oq : 



1 



1/At 



(A2) 



The distribution /i(f2) is normalized so as to yield 



J dSlh(n) = 1. 



It represents the spectral profile associated to the mirror's motion: 8q[£l] = Sqo[h(fl) + h(—Sl)]/2. Therefore, rather 
than the sinusoidal motion of eq. (14), we actually have an exponentially damped sinusoidal motion (where A£ is the 
damping time). In the frequency domain, it amounts to replace the delta function by the distribution h. Since h(fl) 
is very sharply peaked around f^o, however, we may replace ft by Slo every wher e (and thus keep all the development 
of sec. 2), except when computing the vacuum correlation function in eq. ( |Al[ ). For instance, the representation of 
eq. (|q) is replaced by 



5A rc t[a;, u] = exp (ilC(cv)x) 



E 



and moreover fl is replaced by CIq everywhere in «4-i n [0, uj + mfl] except in the argument of the bosonic input operators. 
The finite mechanical linewidth is taken into account only when computing the vacuum correlation function, for which 
we take an average over the distribution h : 



^a in _„(/C_i)a i t nii (/C_i)^ = J d^ J d£l" h(Sl')h(Sl") (Oin| a in „(AC(w + ft'))a in „(AC(w + ft")) 1 |0in) 

The unaveraged correlation function is easily calculated from the commutation relations given by eq. Q 

(0 in| a inn (iC)a inn (/C') f |0 in) = 2ir5(!C - K. 1 ). 



(A3) 



(A4) 



Replacing eq. (A4) into eq. (A3) and using the auxiliary result 

//V 
dnh 2 (n) = — 
2n 

yields 



a ln J/C(^0)«L(^(^- 1 ))) = n^ Z f 



(A5) 



When inserting eq. (A5) into ( Al), it is convenient to employ the variables k x and kn defining a given wavevector 
k and the corresponding bosonic operator a- mn (k x ) [see eq. (^)], rather than the Fourier variable u>. Thus, we take 

k x = —K{lo), and then k = + fc| = — u>. Note that since u> < — k^ eq. ( |l3| ) yields k x > as it should. The 
final result is given by eq. (pfl). 
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APPENDIX B: CONNECTION BETWEEN PHOTON EMISSION RATES AND NUMBER OPERATORS 



In this appendix, we establish the connection between the averaged number operator, as given by eq. ( pq ) and 
the photon emission rate. The field normal mode decomposition as given by eq. (g) is such that the lhs of eq. ( |26| ) 
represents, when multiplied by dk x /{2-n), the average number of (TM polarized) photons with wavevectors whose x 
component is between k x and k x + dk x and whose parallel component is ky . Since the number of allowed values of k|| 
per unit area in the two-dimensional yz reciprocal space is S/(2tt) 2 , the average number of photons with frequency 
between k and k + dk is 

d 3 N = (al tn (k x )a outn (k x ))^d\ L j^. (Bl) 



According to eqs. ( pq ) and (Bl), d 3 N is proportional to the coarse-grained time interval Ai, which allows us to define 
the photon production rate drR = d 3 N/£t. Writing the volume element in spherical coordinates, 

d 3 k = dk x d 2 ku = k 2 dkdn, (B2) 



where Q denotes the solid angle, we obtain from eq. (Bl) the connection between the angular distribution rate of 



photon emission (photon emission rate per frequency and solid angle intervals) and the averaged number operator: 

d 2 R Sk 2 



dkdfl (27r) 3 Af 



a Lt„ (M a out„ {k x )) ■ (B3) 



By replacing eq. ( |26|) into (|B3| ), we obtain the result for the photon emission rate as given by eq. (|28|). 
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Figure Captions 



Figure 1. Photon emission rate as function of the normalized mechanical frequency A (see text), with 6 = 78° 
and kSqo — 0.03. Main plot: as calculated from the perturbation theory (dotted line) and from the nonperturbative 
approach presented in this paper (solid line). Insert plot: numerical results (solid line) and analytical approximation 
where high order sidebands are neglected (dashed line). Logarithmic scale is employed in the vertical axis in the 
insert. Note that for the scale employed in the main plot, it is not possible to distinguish numerical and analytical 
curves. 

Figure 2. Frequency shift A s as calculated from the analytical nonperturbative theory as a function of observation 
angle 0. When 6 varies from to 90°, A s increases by several orders of magnitude. 

Figure 3. wx ku plane representing the field modes. For a given mechanical frequency fio, the region bounded by the 
triangle ABC (dark grey) represents the modes that may be excited. The light grey region represents evanescent field 
modes (not relevant in the present discussion). In this diagram, the photon creation process is associated to frequency 
upshifting. This is represented by a horizontal displacement (since ku is conserved) of length Oq- Accordingly, photons 
in field modes along the BC line (A = 0) are created from vacuum fluctuations that propagate along grazing directions 
(B'C line). 
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